Paramagnetic Spectra of Substituted 
Sapphires — Part I: Ruby* 

By E. 0. SCHULZ-DU BOIS 

(Manuscript received September 29, 1958) 

The paramagnetic resonance properties of Cr ions in AU0 3 (ruby) 
were investigated theoretically and experimentally in order to obtain infor- 
mation necessary for the application of this material as active material in 
a three-level solid-state maser (3LSSM) . Numerically computed energy levels, 
together with their associated eigenvectors, are presented as a function of 
applied magnetic field for various orientations of the magnetic field with 
respect to the crystalline symmetry axis. A more detailed discussion is de- 
voted to energy levels, eigenvectors and transition probabilities at angles 0°, 
54.74° and 90°, where certain simple relations and symmetries hold. Para- 
magnetic spectra for signal frequencies between 5 and 24 kmc are shown; 
agreement between computed and measured resonance fields is satisfactory. 

I. INTRODUCTION 

Among the paramagnetic salts that have been used as active materials 
in three-level solid-state masers (3LSSM), 1, 2 ' 3 ruby shows rather desir- 
able properties. While maser action of this material has been achieved 
at microwave signal frequencies of 3 to 10 kmc, it should be possible to 
cover more than the whole centimeter microwave range. Perhaps even 
more important from a practical point of view are the bulk physical 
properties. Extremely good heat conductivity at low temperatures allows 
handling of relatively high microwave power dissipation. Industrial 
growth of large single crystals by the flame fusion technique and ma- 
chinability with diamond tools make it possible to fabricate long sections 
of ruby to very close tolerances, a necessity in travelling-wave maser 
(TWM) development. Also, ruby can be bonded to metals, thus allowing 
a high degree of versatility in maser structural design. While the use of 
ruby in 3LSSM, in particular in nonreciprocal TWM, will be described 
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in forthcoming papers by members of Bell Telephone Laboratories, this 
paper is intended to give some background on paramagnetic resonance 
behavior of ruby. 

In general, the paramagnetic resonance properties of an ion in a crystal 
can be completely described by a spin Hamiltonian containing a relatively 
small number of constants. In the case of ruby, these include the spec- 
troscopic splitting factors parallel and perpendicular to the crystalline 
axis, ,9 1| and g±, the total spin S = 3/2 and the sign and magnitude of 2D, 
the zero field splitting. Nuclear interactions can be neglected since the 
most abundant isotope, Cr 5z , is nonmagnetic (/ = 0) whereas the mag- 
netic isotope, Cr 53 , (7 = 3/2) has small abundance (9.5 per cent) and 
leads to negligible line broadening only. Taking this into account, one 
can even predict on the basis of the total spin and the crystalline sym- 
metry surrounding the Cr +++ ion that no other terms can occur in the 
spin Hamiltonian. 

However, in order to predict operating conditions of this or other ma- 
terials in a 3LSSM, it is necessary to know the separation of energy 
levels for supplying the proper pump and signal frequencies, the order 
of magnitude of the associated transition probabilities and perhaps other 
circumstances, such as coincidence of transition frequencies, which, by 
spin-spin interaction, may lead to shortening of the associated relaxation 
times (self-doping condition). In this paper, this information is evaluated 
by the formalism of the spin Hamiltonian and, at least in part, compared 
with experiment. The data presented graphically are intended to form 
an "atlas" of the ruby paramagnetic resonance properties. In the paper 
which follows, some general viewpoints are presented on modes in which 
paramagnetic materials can be operated as active materials in a 3LSSM. 
In further papers, paramagnetic spectra of other substitutional ions such 
as Co ++ andFe +++ in sapphire will be presented in order to furnish suffi- 
cient information to find coincidences of transition frequencies of Cr +++ 
with Co ++ or Fe +++ lines resulting in reduced relaxation times (impurity- 
doping condition). 

For a derivation of the method of spin Hamiltonians, reference should 
be made to such review articles as those by Bleaney and Stevens and 
Bowers and Owen. 6 Knowledge of the associated formalism is perhaps 
desirable but not necessary for utilization of the results reported in this 
paper. Briefly, the spin Hamiltonian describes the energy of a para- 
magnetic ion arising from interaction with host crystal environment and 
applied magnetic field. Obeying quantum laws, the ion can exist in one 
of several states associated with discrete energy levels. Transitions be- 
tween such states can occur if the energy balance A/? is supplied to or 
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extracted from the ion. Given some probability for radiative transitions, 
these can be induced by applying a magnetic field of radio frequency 
v = AE/h (h = Planck's constant). If there are more transitions to the 
higher state, net absorption will be observed such as is normally ob- 
served with a spectrometer. If there are more transitions to the lower 
state, stimulated emission of energy will be observed such as is utilized 
for amplification in a 3LSSM. 

II. THE SPIN HAMILTON! AN 

The spin Hamiltonian of Cr +++ in Al 2 :t was first published by Manen- 



kov and Prokhorov , and later by Geusic and Zaripov and Shamonin. 
It was given in the form 

3C = ftHA + gAH x S x + H S ) + D[Si - ±S(S + 1)]. (1) 

The effective spin S = 3/2 is identical with the true spin. All Cr +++ ions 
in the crystal lattice show identical paramagnetic behavior, with the 
magnetic 2-axis being the same as the trigonal symmetry axis of the 
crystal. The best values for the constants seem to be 

2D = -2D' = -0.3831 ± 0.0002 cm -1 = -11.493 ± 0.006 kmc, 

r/,| = 1.9840 ± 0.0000, 

g ± = 1.9807 ± 0.0006. 

While it is customary in spectroscopy to express energy in units of cm 
omitting a factor he (h = Planck's constant, c = velocity of light), units 
of kmc are used simultaneously, omitting a factor of 10 h, because this 
allows direct interpretation in observed spectra. 

In particular, the negative sign of D was obtained by Geusic. He 
deduced this from the fact that g\\ < g ± , since in less than half-filled 
r/-shell ions, such as Cr +++ , the spin-orbit coupling term X is positive, 
and D is given by 2D = \(g\\ — g ± ). Sign and magnitude of D are in 
agreement with results of low-temperature static susceptibility measure- 
ments by Bruger." In this work also, the negative sign of D was con- 
firmed by comparing the relative intensities of two lines at liquid nitrogen 
and helium temperatures. 

The spin Hamiltonian (1) can more conveniently be written in 
spherical coordinates : 

3C = g\\H cos 6S t + \gjl sin 0(c~' v ,S* + + c' v S_) 

- D'\S, 2 - %S(S + 1)]. 
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Here S± = S x =fc iS y . In l>oth representations (1) and (2) the crystalline 
axis was chosen to be the z-axis. While the choice of reference system is 
immaterial to obtaining eigenvalues (energy levels), this choice shows 
up in the associated eigenvectors. The eigenvectors have no direct phys- 
ical interpretation; they must be evaluated in order to obtain transition 
probabilities. The transition probabilities most naturally obtained from 
eigenvectors of the Hamiltonian (2) arc those which correspond to ex- 
citation by rf magnetic fields whose polarization is either linear and 
parallel to, or circular and perpendicular to, the crystalline axis. 

In 3LSSM design, however, it seems more appropriate to analyze the 
performance in terms of rf magnetic fields whose polarization is either 
linear and parallel to, or circular and perpendicular to, the applied field. 
The corresponding eigenvectors and transition probabilities can, of 
course, be obtained from those belonging to the Hamiltonian (2) by a 
4-by-4 transformation matrix. But it is more efficient to obtain them 
directly through a transformation of the original spin Hamiltonian (1) 
or (2) into a coordinate system with the 2-axis parallel to the applied 
field. The result of this transformation is 

3C = (0[ cos 2 6 + g x sin 2 0)j8//& 

- D'(cos 2 6 -\ sin 2 e)[S? - %S(S + 1)] 

- D% cos d sin d[e~ iv (S z S+ + S+S.) + c ,v (&S_ + SJ3,)] (3) 

- D'i sin 2 0(e -2,v S+ 2 + e 2i *S- 2 ). 

III. ENERGY LEVELS AND EIGENVECTORS 

From the Hamiltonian 3C (3), its energy eigenvalues W are found 
numerically by solving the fourth-order secular equation 

II (n I 3C - W | m) || = 0, 

1 (4) 

n,m = 3/2, 1/2, -1/2, -3/2. 

The eigenvalues W are functions of H and 6, but not of <p since, because 
of the symmetry of the Hamiltonian, rotation about the 2-axis does not 
change the physical situation. On the following plots (left-hand sections 
of Figs. 1 through 11) diagrams of energy levels W (in units of kmc) are 
shown as a function of applied field H (in units of kilogauss). Plots are 
given for angles 6 from 0° to 90° in steps of 10° and, in addition, for 
54.74°. 

Also, by change of scales, dimensionless eigenvalues y = W/D' are 
shown as functions of the dimensionless quantity x = G/D', where 

G = (tfn cos 2 6 + g x sin 2 0)&H. 
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Fig. 1 — Energy levels and eigenvectors of the Cr +++ paramagnetic ion in ruby 
Rt Jingle = 0° between crystalline symmetry axis and applied magnetic field. 
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Fig. 2 — Energy levels and eigenvectors at 10°. 
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Fig. 3 — Energy levels and eigenvectors at 20°. 
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Fig. 4 — Energy levels and eigenvectors at 30°. 
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Fig. 5 — Energy levels and eigenvectors at 40°. 



FIELD PARAMETER, G/D 
12 3 4 

1 




2 4 6 8 10 

APPLIED MAGNETIC FIELD, H, 
IN KILO-OERSTEDS 



2 3 

FIELD PARAMETER. 



G/D 



Fig. 6 — Energy levels and eigenvectors at 50°. 
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Fig. 7 — Energy levels and eigenvectors at 54.7°. 
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Fig. 8 — Energy levels and eigenvectors at 60°. 
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Fig. 9 — Energy levels and eigenvectors at 70°. 
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Fig. 10 — Energy levels and eigenvectors at 80°. 
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Fig. 11 — Energy levels and eigenvectors at 90°. 



This dimensionless representation facilitates computations and reveals 
more clearly symmetries and singular relations in the energy level 
scheme. It also permits the use of the same diagrams for ions having the 
same Hamiltonian but different zero field splitting 2D. Similar energy 
level diagrams were computed by P. M. Parker 11 for the case of nuclear 
spin resonance with nuclear quadrupole splitting present which is de- 
scribed by the same type of Hamiltonian. 

As a convenient way to identify the energy levels IF, a quantum num- 
ber n ranging from — f to +f is used in order of increasing energy. 
Thus W(— f) is the lowest, TF(f) the highest energy level. It is easily 
shown that, for all angles 8 and x = 0, y{ — J) = y(—?) = — 1 and y(?) 
— 2/(1 ) = 1 • As a matter of mathematical curiosity, it may be mentioned 
that, irrespective of 8 at x = 1, y(?) = 1/2. 

The eigenstates | ft) (using Dirac's "ket" notation) associated with 
energy levels W(n) can be expanded in the form 



3/2 

»> = 2 a(n; m) I m). 

m— 3/2 



(5) 



Here, | m) are eigenstates of a Zeeman Hamiltonian 3C = gfiHS, . The 
x(n; m) are amplitudes of eigenvector components or, more briefly, 
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eigenvectors and form a normalized and orthogonal system of coefficients. 
With high applied magnetic field II, | n) — * | n) and a(n; n) — » 1; there- 
fore, | m) are termed high-field eigenstates and «(fl; m) high-field eigen- 
vectors. 

Eigenvectors a{n; m) are obtained as solutions of linear homogeneous 
equation systems, the matrix of which forms the secular equation (4) 
with the particular eigenvalue W(n) inserted. Since this matrix depends 
on tp, the a(n; m) are also functions of tp. The computations were carried 
out for tp = 0, with 9 restricted to < 6 < tt/2 and negative sign of 
D = -D'. 

This choice implies that the crystalline axis lies in the positive quad- 
rant of the x-z plane and it results in real eigenvectors a(n; m). These 
are plotted in the right-hand sections of Figs. 1 through 1 1 , adjacent to 
plots of the corresponding eigenvalues W(n). Negative a(n; m) are indi- 
cated by dashed lines. 

A nonzero <p would, in general, result in new complex eigenvectors 
a'(n; m) = [exp i(m - n)<p]a(n; m). Taking tt/2 < d < t or tp = ir 
would change the sign of every second eigenvector, that is, of those with 
m as n ± ] and m = n ± 3. The same is true for a change of sign of D, 
but then, in addition, every n and m and energy eigenvalue has to be 
replaced by its negative. It is obvious that such transformations do not 
change the physical situation as far as transition probabilities are con- 
cerned. 

IV. TRANSITION PROBABILITIES 

There are several ways in which transition probabilities could be 
evaluated and plotted. < >ne way would be to consider transitions induced 
by radiation of given polarization. With eigenvectors belonging to the 
Hamiltonian (3), the obvious rf magnetic field polarizations to consider 
are those with rf //-field linear and parallel to, or circular and per- 
pendicular to, the applied field. But transitions due to any other polariza- 
tion could be evaluated as well. Perhaps more natural from a theoretical 
point of view would be an evaluation of the maximum transition prob- 
ability. This requires a particulai — in general — elliptical, polarization for 
excitation, which of course should be evaluated, too. All polarizations 
orthogonal to this (which in general are elliptical as well), and which 
describe a plane in space having complex components, are associated 
with zero transition probability. Taking into account these different 
viewpoints and the six transitions which are possible between four 
energy levels, it appears that an unrealistically high number of graphs 
would be necessary to describe the transition probabilities properly. 
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Furthermore, in maser design it is usually sufficient to know the order 
of magnitude of transition probabilities of particular lines, because often 
other factors may be more important. Therefore, no plots of transition 
probabilities are presented. On the other hand, enough of the per- 
tinent formalism is given below so that any transition probability 
can be evaluated from the eigenvectors plotted. 

Following essentially Bloembcrgen, Purcell and Pound, 1 "' with slight 
generalization, the transition probability w describing the rate of transi- 
tions per ion from a lower state n to a higher state n' > /7 is given by 



ir 



M - = i(^J *(" - "»)l <n\&\fi)\*. 



(6) 



Here Hi is the amplitude of the exciting rf magnetic field, g{i> — v ) is a 
normalized function describing the line shape fg(v — vq) dv = 1, and S\ 
is a spin operator reflecting the polarization of the inducing rf magnetic 
field. If the rf magnetic field is described by the real parts of 
H x = Hide*" 1 , H„ = Hibc' ", H z = Hice"" with "complex direction 
cosines" a, b, c accounting for elliptical polarization, 

a*a + b*b + c*c = 1, 

then 

S } = a*S x + b*S u + c*S z . (7) 

Matrix elements for »S'i occurring squared in (G) are linear combinations 
of the following three: 

(n | S z | ft) = 2 ma{n; m)a(n; m), (8) 

m— 8/2 

(n I 8 + | ft) = 

+1/2 (9) 

2 [S(S + 1) - (m + l)m] u2 a(ii'; m + l)a(/7; m), 
m— 3/2 

(n | SL | ft) = 

+3/2 (10) 

£ [S(S + 1) - (m - l)m] m a(n; m - l)a(/7; m). 

m— 1/2 

The square root in (9) and (10) takes on the values \/s, 2 and \/3. 
For example, with linear polarization in the ^-direction, H : = Hi cos w/ 
and S t = S z . For circular polarization perpendicular to the ^-direction. 
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//, = (l/y/2)Hi cos at, H y = ±(l/\/2Wisinco/and & = (l/V2)S±. 
For linear polarization in the x direction, H z = H\ cos bit and Si = S x = 
%(S+ + S-). Similarly, in the y direction H„ = H x cos bit and & = 

,s, ; = (i/ao(s+ + «-). 

The expression (7), or more correctly, the associated matrix element, 
can be interpreted as a scalar product of (a*, b*, c*) with (ft! \ S \ ft). 
It should he noted that, in general, all components can be complex. As 
a consequence of this interpretation, the maximum transition probability 
occurs if II r! or (a, b, c) is parallel in space and conjugate complex in 
phase to (ft' \ S \ n). Since for real eigenvectors the matrices (8), (9), (10) 
are all real, it follows that (ft! S x \ n) and (n' \ S 2 \ n) are real, whereas 
<n' | S y j ft) is imaginary. Thus, for all ruby lines, the polarization for 
maximum transition probability will be a linear combination of H x and 
//, components with an H y component in quadrature. In a similar fash- 
ion, a set of complex direction cosines can be found which causes the 
scalar product of (a*, b*, c*) with (ft' \ S\n), and hence the transi- 
tion probability, to vanish. These vectors (a, 6, c) describe a plane orthog- 
onal to the vector for maximum transition probability. 

It should be noted that frequently the complete formula (0) is not used 
to evaluate and compare transition probabilities. Instead, usually only 
the squared matrix element | (ft' \ S t \ ft) | 2 is computed and this is then 
compared with a simple standard transition. The obvious standard is 
the transition -1/2 -* +1/2 of an S = 1/2 Zeeman doublet induced 
by circular polarization. This is described, in our notation, by 
| ( -fi/2 | (1/V / 2)'S' + | —1/2) | 2 = 1/2. Accordingly, transitions involv- 
ing a squared matrix element of order 1 or greater are considered strong, 
while perhaps 1/100 is typical of weak transitions. 

V. SPECIAL CASES 

5.i. e = o°. 

The energy levels are parts of straight lines y = 1 ± hx, — 1 ± \x 
with change of slope for some of them at x = 1 and 2. Eigenvectors 
are ±1 and only, again joined for some levels at x = 1 and 2. The 
minus sign of eigenvectors at 0° has no significance; it is only used to 
preserve continuity to neighboring angles. 

At 6 = 0° and .r < 2, the labeling of energy levels by high field quan- 
tum numbers in order of increasing energy is perhaps not the usual one. 
In this paper, lftnvever, it seems appropriate because, with this termi- 
nology, in going from 6 = 0° toother orientations, the notation of states 
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stays the same. It may be pointed out that energy levels defined in this 
fashion should be considered as continuous functions of applied field 
without cross-overs (see Fig. 1). The reason is that any off-diagonal 
perturbation will indeed prevent levels from intercepting by perturbation 
theory arguments. 
Only three transitions are allowed: 

< x < 1 : (+1 1 S + | +J> 2 = 4 

<+I I S- | -|) 2 = <+| | S + | -I) 2 = 3, 

1 <x< 2: <+!|£ + | -I) 2 = 4 

(+1 1 s- 1 +W = <+i I s + \ -|) 2 = 3, 

2 <x: (+I|<$'+I -|> 2 = 4 

(+1 I S + | +|> 2 = <-?|5+|-!> 2 = 3. 

It is interesting to note that, for < x < 2, one transition requires 
opposite polarization from the others. This was verified in an experi- 
ment. Resonance absorption was measured for this and another transi- 
tion in a propagating comb-type slow-wave structure having regions 
of predominantly circular polarization. Reversal of applied magnetic 
field results in drastic increase of one and reduction of the other line. 

5.2. e = 54.74°, cos 2 e = i/s. 

For this angle, the fourth-order secular equation reduces to a bi- 
quadratic one. The four eigenvalues are y = ±[1 + f* ± (3a; 2 + 
x) "]' 2 . This implies an up-down symmetry y( — ii) = —y{n). The 
closest approach of the two middle eigenvalues is y(+i-) — y{— i) = 
1 t at x = 1. A similar symmetry relation holds for eigenvectors 
a{ — n; —m) = (nm/\ nm\)a(n; in). Asa consequence, some transition 
probabilities for linear polarization are identical, namely 

<-J I S. | -I) = <+! | 8, 1 +|) 
and 

(+5I&I-5)- -<+5l&l-t>. 

The analogous is not true for other polarizations. 
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5.3. 9 = 90°. 

The secular equation can be factorized into two quadratic equations 
with the solutions 



2\l/2 



2/(1) = % + (1 + * + xV\ 



y(^)= -|+(i-* + *T 1 , 



iK-i) - 5 - (i + * + *r J , 



*(-!) = -|-(i-* + * 2 ) 1/2 . 

Each state contains only two eigenvectors, namely a(n; n) and 
a{n; n ± 2). In addition, a(w; n) = a(n ± 2;w ± 2) and a(n; n ± 2) = 
— a(n ± 2; w). As a result, transition probabilities between adjacent 
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Fig. 12 — Paramagnetic resonance spectrum of Cr +++ ions in ruby at signal fre- 
quency 5.18 kmc. 
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Fig. 13 — Resonance spectrum at 6. OS kmc. 
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ANGLE, 6, WITH RESPECT TO CRYSTALLINE AXIS 

Fig. 14 — Resonance spectrum at 9.30 kmc. 
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levels n —* n + 1 contain only matrix elements of S + and <S_ , the same 
being true for -| -» +|. Double jumps w -» w + 2 are described by 
nonvanishing elements of S, only. 

VI. PARAMAGNETIC RESONANCE SPECTRA 

In Figs. 12 through 17 some resonance spectra are shown for signal 
frequencies of 5.18, 6.08, 9.30, 12.33, 18.2 and 23.9 kmc. The plots show 
resonance fields as functions of the angle between crystalline axis and 
applied field. Measurements have been carried out at all of these fre- 
quencies to varying extents, although measured values are recorded only 
on Figs. 14 and 15. Generally, these spectra have been used in the lab- 
oratory to align ruby crystals by resonance for maser experiments. They 
have proved accurate to about ±50 gauss. 
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ANGLE, 6, WITH RESPECT TO CRYSTALLINE AXIS 

Fig. 15 — Resonance spectrum at 12.33 kmc. 
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ANGLE, 9, WITH RESPECT TO CRYSTALLINE AXIS 



Fig. 16 — Resonance spectrum at 18.2 kmc. 

Measurements at 9.3 kmc are an extension of Geusic's work 4 and con- 
firm his results. Results at 12.33 kmc show some discrepancy between 
theory and experiment, which, however, is believed to be caused by 
inadequate magnetic field measuring equipment used in an experiment 
designed for other purposes. As a general rule, the spectra show two 
looping lines if v < 2D. Lines marked "forbidden" are strictly forbidden 
at 0° only. Usually, however, they can be followed quite close to 0° by 
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Fig. 17 — Resonance spectrum at 23.9 knic. 
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use of more sensitivity in the spectrometer. An exception is the line 
shown on the graphs having the lowest resonance field at 0° if v < f D'. 
It has the second lowest resonance field if f D' < v < f /)' and the third 
lowest if f // < v < 3D'. It originates between — f and +$ eigenstates 
at 0° and is more strongly forbidden than the other forbidden lines; hence 
it usually ceases to be measurable at about 30°. 

For reasons of symmetry, all lines approach 0° and 90° with zero slope 
dH /dd. Experimentally, it has been found that most lines are rather 
narrow at 90° and similarly at 0°, whereas they broaden in proportion 
with dH/dd. This behavior is expected from crystalline imperfections if 
these can be interpreted as fluctuations throughout the crystal of the 
direction of the crystalline axis. 
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